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Abstract 

A class of recently introduced su{n) 'free-fermion' models has recently been used 
to construct generalized Hubbard models. I derive an algebra defining the 'free- 
fermion' models and give new classes of solutions. I then introduce a conjugation 
matrix and give a new and simple proof of the corresponding decorated Yang-Baxter 
equation. This provides the algebraic tools required to couple in an integrable way 
two copies of free-fermion models. Complete integrability of the resulting Hubbard- 
like models is shown by exhibiting their L and R matrices. Local symmetries of the 
models are discussed. The diagonalization of the free-fermion models is carried out 
using the algebraic Bethe Ansatz. 
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1 Introduction 



The two-dimensional Hubbard model was introduced to describe the effects of correlation 
for d-electrons in transition metals It was then shown to be relevant to the study of 
high-Tc superconductivity of cuprate compounds. 

The one-dimensional version also has interesting features. The model is integrable 0, 
^, ^ . The integrability framework of the model is the quantum inverse scattering method 
P, ^ 1^. However, despite sharing many properties with discrete quantum integrable 
models, the model had a peculiar integrable structure which defined a class of its own. It 
was therefore natural to look for integrable generalizations. Mapping the fermionic model 
to a bosonic one with a Jordan- Wigner transformation reveals interesting properties. The 
local fermionic symmetries become non-local. This has been known for some time but 
does not seem to have been explored further [^. 

Another interesting feature of the ID Hubbard model and of most interacting one- 
dimensional systems is their Luttinger liquid behavior Such a behavior is not 
restricted to one dimension however ||10|| . 

A bosonic su{n) Hubbard model which contains the usual su{2) model was recently 
introduced in |T1[]. These models were shown to be integrable and to have an extended 



su{n) symmetry [|I2|. The model is built by coupling two copies of the recently discovered 
su{n) XX 'free-fermions' models For n = 2 a fermionic formulation exists, but for 

n > 2 finding an analogous framework remains a tantalizing problem. Strictly speaking, 
the name free-fermions model does not seem appropriate. It is easy to convince oneself on 
dimensional grounds that an expression in terms of fermionic operators can only happen 
for a subclass of models. However for want of a better characterization I shall stick to the 
foregoing appellation. 

Other types of generalizations of the Hubbard model exist. They are mostly of the 
fermionic type, that is, built from fermionic operators; see for example [0, [T^, . 

In this work I look for new solutions of the Yang-Baxter equation which share the 
same features as the known XX models. I derive an algebra which unifies the different 
'free-fermions' representations, and greatly simplifies the calculations. This algebra is 
reminiscent of the Temperley-Lieb algebra. The former algebra is more restrictive and all 
the representations found so far are also representations of the Temperley-Lieb algebra. 
Defining a conjugation operator allows for a simple and new derivation of the decorated 
Yang-Baxter equation. This equation, introduced by Shastry while studying the usual 
Hubbard model, is an important algebraic component of the integrability proof for the 
'fusion' of two models. The 'fusion' or coupling of two commuting free-fermions copies is 
then described, along with an algebraic proof of the integrability of the resulting Hubbard- 
like models. I then give new representations of the free-fermions algebra and of the 
conjugation matrix. I discuss symmetry issues related to these models. Diagonalization 
of the 'free-fermions' models using the algebraic Bethe Ansatz method, shows that their 
spectrum is highly degenerate and simple in a certain sense. Some outstanding issues and 
possible directions are discussed in the conclusion. 



1 



2 A new algebra 



Let -E"^ be the n x n matrix with a one at row a and column and zeros otherwise. 
Consider the i?-matrix of the su{n) XX model ||13| : 



R{\) = a(A) [E"" O E"" + ® ^^^] 

ci,/3<n 



+c(A) ^ ® E'^" + E"" O E"") 



a<n 



where a(A) = cos(A), b = sin(A) and c(A) = 1. The functions a, b and c satisfy the 
'free-fermion' condition: + 6^ = c^. For this set of parameters, a Jordan- Wigner 
transformation turns the U = hamiltonian density for su{2) into a fermionic expression 
for free fermions hopping on the lattice. 

I now look for i?-matrices having the above form, namely 

R{X) = p(i) + p(2) cos(A) + P(3) sin(A) (2) 
and impose the property of regularity 

^(0) = I (3) 

where I is the identity operator. One therefore has P^^-* + P'-^^ = I. There is no loss 
of generality in choosing the proportionality constant to be one since a solution to the 
Yang-Baxter equation is defined up to a multiplicative function of A. 
Requiring R to satisfy the Yang-Baxter equation 

i?12(A)P23(A + fi)RM = P23(/^)^12(A + /X)P23(A) (4) 

yields a finite set of equations. One develops on a set of hnearly independent functions 
and equates the operatorial coefficients. I do not reproduce here all the equations and 
rather concentrate on the basic ones: 

p(2) + (p(3))2^^j (5) 

[A + B, [A, B]] = (1 - 3a){A - B) + A^ - B^ (6) 
A^ + {1- 2a)A^ = (3 I (7) 

(3) (3) 

where A = P12 , B = P23 and a, (3 are two arbitrary complex numbers that arise upon 
solving equations of the type M12 = M23. Equation and the regularity equation allow 
to keep P'-^-* as the sole unknown operator. All other equations are therefore equations 
for P^^\ Equation (BI) means that the Reshetikhin criterion is satisfied ITS, IT3. This can 



be seen as an integrability test for spin chains pO[ . It seems at first that the large number 
of constraints on one operator cannot have a solution. But we already know that (|l]) is 
a solution for which a = 1 and (3 = 0. I have looked for the minimal set of equations 
which solves all the equations involved in the Yang-Baxter equation and have found the 
following algebra. 
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(3) 

Let Ei = Pjj^i, that is, Ei acts non-trivially on the adjacent spaces i, i + 1. The 
defining relations of the free-fermions algebra A are: 

{Ef , Ei±i} = Ei±i , Ef = Ei (8) 
EiEi^.Ei = , EiEj = EjEi for \t-j\>2 (9) 

where {A, B} = AB + BA. The fourth equation just expresses the fact that E^ and 
Ej commute when they act non-trivially on disjoint spaces. There does not seem to be 
solutions of the Yang-Baxter equations for the foregoing i?-Ansatz unless a = 1 and 
/3 = 0. 

The above algebra is reminiscent of the Temperley-Lieb (TL) algebra All the 

solutions given in section ^ can be put, after a 'gauge transformation' (a kind of simi- 
larity transformation), in a Temperley-Lieb form. Conversely, the set of solutions of the 
Temperley-Lieb algebra is much larger and most of its solutions do not map to a free- 
fermions system. Thus the algebra A is much more restrictive. More details are given in 
section H. 

One then checks that P'-^^ and P'-^^ form a complete set of projectors on the tensor 
product space C" ® C" 



p(l) _^ p(2) ^ J ^ (p(l))2 = p(l) ^ (p(2))2 ^ p(2) ^ p(l) p(2) ^ p(2) p(l) ^ q ^^q^ 

The operator P^^^ is a square root of the operator P'-^\ and p(2)p(3) = p(3)p(2) = g. 
The above relations imply that the matrix R satisfies the unitarity property 

P(A)P(-A) = I cos^ A (11) 



3 The decorated Yang-Baxter equation 

The decorated Yang-Baxter equation is an equation similar in form to the Yang-Baxter 
equation. It was first introduced in as an algebraic relation at the root of the inte- 
grability of the su{2) bosonic Hubbard model. This underwent a first generalization in 
|1 1[ | . I give here the ingredients needed for the existence of the DYBE. 

Assume there exists a 'conjugation' matrix C acting on one copy C" such that 

c^ = i, {a,p;f} = o , z = i, 2 (12) 

CiPi? = Pi?C^2 , 2(pf )2 = I - (13) 

where Ci = C®I and C2 = I® C. Then Ci commutes with P^^2^ and Pil"*- These relations 
imply the following, equivalent, conjugation relations for R and R = VR, where V is the 
permutation operator on C" C": 

CiRuiX) = Pi2(-A)a , ^ = 1, 2 (14) 
QRuiX) = Pi2(-A)C, , 2, J = 1,2, t^j (15) 

Note that only the anticommutator relations were used here. The remaining relations will 
be needed in section [4.2| . For the P-matrix (|l]) one has C = J2a<n ~ E^^. 
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That the i?-matrices satisfy a decorated Yang-Baxter equations (DYBE) is a simple 
consequence of the conjugation relations and of the Yang-Baxter equation. Consider the 
following version of the latter equation: 

RuiX - /i)i?13(A)i?23(/i) = RMR23iX)Rl2iX - fl) (16) 

One multiplies ([T6|) by C1C2, commutes them appropriately and removes one C using 
= I.Q After letting fi — > — /i one obtains the following version of the DYBE 

Ri2{X + /i) Ci i?i3(A) R2M = RM i?23(A) C2 Ri2{X + fi) (17) 

It is worth noting that, while the YBE is invariant under a 'gauge transformation', the 
DYBE is not. This is due to the arguments of the matrices. 

At this point, it is not clear whether the existence of C follows from that of P^^\ 
although it is the case for the models of section ^. 

In [ITTl , |T^ we saw how to couple two su{n) XX models a la Shastry. We generalize 
this procedure to the foregoing models. 



4 A new kind of fusion 

Coupling two su{2) XX models in an integrable way in the framework of the Quantum 
Inverse Scattering Method was first done by Shastry It was then generalized in 
| TT[| . The results of the above sections provide all the ingredients required to couple 
two solutions of the algebra A. Thus the derivation of the two following subsections is 
algebraic and does not depend on a specific representation. 

In the context of the quantum inverse scattering method and i?-matrices, the word 
fusion has a precise meaning. It refers to the construction of higher-dimensional solutions 
to the Yang-Baxter equation using lower-dimensional solutions. One uses projection op- 
erators and reduces the resulting fused-space dimensionality. For instance, fusing two 
spin-^ su{2) i?-matrices results, after projections, in an i?-matrix for the spin-lxspin-1 
representation, with dimension 3^ x 3^. 

In the kind of fusion described below, one couples two models without a reduction in 
dimension for the tensor space of the new model; no projection is implemented. The term 
fusion here then takes a different meaning. 



4.1 Lax and transfer matrices 

The transfer matrix is the generating functional of the infinite set of conserved quantities. 



The construction given in [jTT| is still valid here. We consider two commuting copies of the 
free-fermion models found in the preceding section. Let us stress that the unprimed and 
primed copies need not be of the same type. For instance, the 'left' copy can be (ni, 77-2) 
while the 'right' copy is (?t-i,?t-2) with n not necessarily equal to n (see section ^ for the 
definitions and notation). The Hilbert space of the chain is a tensor product of C" *^ , the 
Hilbert space at each site. 

^This simple derivation of the DYBE does not appear in and appears to have been overlooked in 
the literature. 
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Consider the matrix 

Io{h) = cosh(|) I + sinh(|) Co C; = exp (^^ Co C',^ (18) 

where C and C' are the conjugation matrices of the models being coupled. The second 
equality follows from = I. The parameter h is related to the spectral parameter A by 

sinh(2/i) = [/sin(2A) (19) 

where U characterizes the strength of the coupling. One chooses for h{X) the principal 



branch, which vanishes for vanishing A or U. Then for U = the monodromy matrix (21) 
becomes a tensor product of the two, uncoupled, models. The Lax operator at site i is 
given by: 

Lo.(A) = Io{h) Ro^{X) i?;,(A) Io{h) (20) 
and the monodromy matrix is a product of Lax operators 

T(A) = Lom(A)...Loi(A) (21) 

where M is the number of sites on the chain. The transfer matrix is the trace of the 
monodromy matrix over the auxiliary space 0: 

r(A) = Tro [(Lom-^oi) (A)] (22) 

The conserved quantities are then given by 

/rfPlnr(A)\ , , 

Hp+i = I 1 ' P > (23) 



dX'P 



A=0 



This completes the 'fusion' of the two models. 

A proof that H2 commutes with r(A) was given in [0. This proof is algebraic and 
holds also for the models considered here. It yields relation (0). The proof in the 
following section ensures the complete integrability and yields ([T9| ) again. 

The derivative of the matrix / gives the coupling term appearing in (|45|). Note that the 
definition involving a logarithm has two benefits. Besides giving the most local operators, 
it further disentangles the two copies. 

The construction of a non-additive i?-matrix intertwining two L-matrices at different 
spectral parameters goes through exactly as in . This automatically implies the exact 



integrability of the models with periodic boundary conditions. 
4.2 The i?- matrix 

The i?-matrix which intertwines the Lax operators, 

i?(Ai, A2) L(Ai) ® L(A2) = L{X2) ® L(Ai) i?(Ai, A2) , (24) 

is again given by 

i?(Ai,A2) = /l2(M/34(/^l)(«i?13(Al-A2)i?24(Al-A2) 

+/3^13(Ai + A2)Ci^24(Al + A2)C2) Il2{-hl)l34{-h2) (25) 
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This matrix acts on the product of four auxihary spaces labeled from 1 to 4, and a, P are 
to be determined. One then requires relation (0) to be satisfied and uses (p!6|) and ( |T7|) 
to derive the following equation: 

(«/?i3(Ai - A2)/?24(Ai - \2) + P CsRniXi + \2)C^R2^{\i + As)) /12 (2/11)134(2/12) = 

/l2(2/l2)/34(2/ll) (ai?i3(Ai - A2)/224(Al - A2) +/3i?l3(Ai + A2)Cii?24(Al + A2)C2) 

Expanding the exponentials (see (p!8|)) and using (0) results in the cancellation of half 
the terms on each side. Using (|^) and relations (p!2|-[l3D for all the terms, yield only two 
equations: 

— = — tanh(/ii + /i2) , — = — tanh(/ii — /i2) (26) 
a B a A 

where a = cos(Ai — A2), b = sin(Ai — A2), A = cos(Ai + A2) and B = sin(Ai + A2). The 
compatibility equation 

tan(Ai - A2) _ tanh(/ii - /is) 
tan(Ai + A2) tanh(/ii + /i2) 

is satisfied if and only if equation (p!9D is satisfied for the pairs (Ai, hi) and (A2, h2)- One 
can then factor out a = a(Ai,A2) which appears as an arbitrary normalization of the 
i?-matrix, to obtain: 

^(Ai, A2) = a(Ai, A2)/i2(/i2)/34(/ii) ( i?i3(Ai - A2)i?24(Ai - A2) + "'""i^ 7 M 

y sm(^Ai + A2) 

X tanh(/li + h2)RlsiXl + A2)Cii?24(Al + X2)C2) /i2(-/ii)/34(-/^2) (28) 



The underlying algebraic structure at work here is the one elucidated in sections ^ 
and 0. We stress here that this proof is rigorous and valid independently of the specific 
representation for P^^^ and C. It only involves algebraic properties. 

We conclude the integrability proof. The monodromy matrix being a tensor product 
of M copies of L matrices, one has 

i?(Ai, A2) T(Ai) ® T(A2) = T(A2) ® T(Ai) i?(Ai, A2) (29) 

Taking the trace over the auxiliary spaces, and using the cyclicity property of the trace, 
one obtains [r(Ai), r(A2)] = 0. We have thus proven that all the conserved charges Hp 
mutually commute. 

The matrix R satisfies the regularity property 

^(Ai,Ai) = a(Ai,Ai) I (30) 

and the unitarity property: 

^(Ai,A2)i?(A2,Ai) = a2(Ai,A2)cos2(Ai - A2) 

X (^cos^(Ai - A2) - cos^(Ai + A2) tanh^(/ii - /i2)) I (31) 

The intertwiner R satisfies a Yang-Baxter relation of its own: 

i?12(A2, A3)i?23(Al, A3)i?i2(Ai, A2) = -R23(Al, A2)-Rl2(Al, A3)i?23(A2, A3) (32) 
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where A and h are related through (|T^). As explained in |]T2| the direct verification of 
this relation is tedious, but can be avoided. The proof done for su{2) in generalizes; 
this proof is based on Korepanov's tetrahedral Zamolodchikov algebra EM. All variants 
of the proof follow from the foregoing algebraic structure. A notable feature of the matrix 
-R(Ai, A2) is its non-additivity property; the A dependence cannot be reduced to a difference 
(Ai — A2). This is the source of the difficulty in verifying p^). 

I now give new solutions of the YBE which have the foregoing properties. 



5 New models 

Let n, rii and ^2 be three integers such that 

2<n, l<?T-i<'n,2<?T. — 1, rti + n2 = n (33) 

and A, B be two disjoint sets whose union is the set of basis states of C", with card(y4) = rii 
and card(i?) = n2. Let 

p(3) = i^a^E'^" ® + ^a^E"^ ® E^") (34) 

a6A/3g_B 

p(i) = (p{3))2 =J2J2 {e^" ® + ® E^^) (35) 

P(2) _ jj _ p{l) _ ^ ^'^^ ® E'*''*' + ^ E/3'/5' (36) 

a,a'eA I3,P'&B 

The ni.n2 parameters Xap are arbitrary complex numbers. Latin indices belong to A while 
greek indices belong to B. These operators satisfy all the defining relations of the algebra 
A and therefore one has an i?-matrix 

R[\) = pW + p(2) cos(A) + P^^^ sin(A) (37) 

which satisfies the YBE. Denote this representation by (ni,n2). 

The conjugation matrix is defined up to an overall sign; it is given by 

C =Y,E^^ -Y, E^^ (38) 

and satisfies all the relations of the preceding section. The DYBE therefore holds. 

Note that, unless one wants to perform numerical calculations or write down an x 
matrix representation, one need not specify which states belong to which set, thereby 
keeping a 'symmetrical' labeling. Note also that the restriction rii < n2 is not essential. 
It just avoids a double counting of distinct models since one has the obvious symmetry 
A ^ B. The number of models, for a given n, is equal to the integer part of n/2: [|]. 

It is possible to perform a gauge transformation which puts these models in a TL 
form.0 One obtains 

if^^X) = I cos(A) + (P(=^) + pTL) sin(A) (39) 
pTL _ iY^Y. {e^^ ® - E"" ® Ef^^) (40) 

a£Al3eB 

^The connection between some XX models and their TL formulation was pointed out by Martins p5{. 



7 



where i'^ = —1. The operator E = P^^^ + P^^ satisfies the Temperley-Lieb algebra 

E^ = 0, E,E,±iEi = E, , EiEj = EjE, for \t-j\>2 (41) 

and R^^{X) the Yang-Baxter equation. Let M be any invertible matrix. A general class 
of solution of the TL algebra is given by 

Eab,cd = MabM-/ , tr(*M.M-i) = (42) 

where, on the left-hand side, a, b (c, d) are the row (column) indices. However, only for 
n = 2 and M oc antidiag(l, ±z) does one have a free-fermions model. The other {ni^n2) 
models do not fit in this scheme. This confirms the statement made in section ^ 

For ni = 1, n2 = n — 1, and all the Xap equal to each other, one recovers the su{n) 



XX models found in ||T3[. Allowing the twist parameters x to be unequal amounts to 
a multiple deformation of these models. We see in the next section that the degree of 
symmetry depends on the x's. 

Before diagonalizing the conserved quantities of the free-fermion models we pause to 
consider their symmetries and their quadratic defining hamiltonians. 



6 Hamiltonians and symmetries 

Periodic boundary conditions are assumed. Consider first the free-fermions models. The 
quadratic hamiltonian calculated by ([23| ) is given by 

^2 = E 4?i (43) 

i 

The cubic hamiltonian is equal to 

^3 = E[45i,^£U] (44) 



This is a boosted form of H2 ||2^. The commutation of H2 and is a consequence of 
the Reshetikhin criterion (||). We expect considerations about the boost structure and 
the explicit form of the conserved quantities to generalize straightforwardly |1^ . 



The quadratic Hubbard-like hamiltonians obtained by fusion are given by 

H2 = Y.pSIi + T.P^i + UT.C^^ (45) 

i i i 

where primed and unprimed quantities correspond to the two commuting copies. The 



cubic hamiltonian is not given by a boosted form of H2 [[Til , p7 |. 

The hamiltonians H2 are defined in one dimension but can be evidently defined on 
any lattice; integrability is lost however. These hamiltonians can be written simply in 
terms of su{n) hermitian traceless matrices. For \xai3\ = 1 and U real the hamiltonians 
are hermitian. 

Because of the structure of the hamiltonians one expects, at least, the diagonal gen- 
erators to commute with all the conserved quantities. As seen for the models of |1^ the 
symmetry may be larger. One has the following relations 

Va,6eA [^f,Lo.(A)] = -[Eo"^Lo.(A)] m x^f, = x,p ^ (3 e B (46) 
Va,/3G5 [Ef,Lo^{X)] = -[Ef,U,{\)] iff x.^ = x.^j Va G A (47) 
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where L = R. The hnear magnetic-field operators 



= Y.Ef a.heA (48) 

i 

Hf = Y.Ef a,(3eB (49) 

i 

commute with the transfer matrix if the parameters Xap satisfy the above constraints for 
the corresponding indices. One just uses the expression (|2^) , its cychc structure and the 
relations (|6) and (^Tj). The commutation with r(A) implies the commutation with all 



the hamiltonians Hp. 

In particular all the diagonal operators are symmetries, without any constraints on the 
parameters x. When all the parameters are equal to one parameter, say Xap = x , Va G 
A , V/9 G B, the full local symmetry is su{ni) © su{n2) © w(l). This symmetry is largest 
(in terms of number of generators) for ni = 1 and n2 = n — 1, that is for the models of 
reference . 



It is straightforward to conclude that the models obtained by fusion inherit the lo- 
cal symmetries of their components. Again, when — ^^'^ ^^e left copy and 
~ f*^^ right copy, the full local symmetry of the model (ni,n2) x {n^^n^) 
is {su{ni) © su{n2) © m(1)) x {su{n^) © su{n2) © m(1)) . 

One consequence is that one can add magnetic-field terms for each symmetry genera- 
tors, without spoiling the integrability of the models. 



7 Algebraic Bethe Ansatz 

The diagonalization by Bethe Ansatz of the free-fermions models is very similar to the one 



for the su{n) XX models. See [13] for additional details omitted here. I am considering 
the to avoid unessential complications. 

The transfer matrix defined in (0) is the trace over the auxiliary space of the mon- 
odromy matrix T(A). The latter is an n-dimensional matrix whose entries are operators 
acting on the Hilbert space C" © ... © C", with a copy for every site. The number of 
sites is M. Some elements of the monodromy matrix are used to create Ansatze for the 
eigenstates. When written in components, equation (^) provides the algebraic relations 
needed to find the action of the transfer matrix on the states. 

We now use the following notation for some elements of the monodromy matrix: 

S = Tii, Ca = Tia , a = 2, = Ti^ , p = ni + l,...,n (50) 

where, as usual, the Latin indices belong to A and the greek indices to B. The remaining 
elements are denoted by T^,^. The transfer matrix is given by 

rii n 

ri\) = Si\) + Y.Taa{\)+ E TppW (51) 

a=2 /3=ni+l 

It is easy to see that the vector ||1) = |1) © ... © |1) is an eigenvector of the transfer 
matrix. The only non- vanishing elements of the monodromy matrix on ||1) are: 

Six) 111) = (cos(A))"'^||l) , TppiX) 111) = (x-Sin(A))"'^||l) (52) 
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and the action of all the C operators, Ca (a 7^ 1) and C/3. 

It turns out that it is still possible to construct Bethe Ansatz eigenvectors using the 
C/3 only, namely: 

|Ai,...,Ap)= J2 i^"--'"^^^ai(Ai)...C,,(Ap) 111) (53) 

«p,...,ai 

where the parameters Aj and the coefficients F are to be determined. This Ansatz also 
vanishes identically if p > M; the proof of this fact is given below. Equation (^9]) gives 
the following relations: 

CMCM = C,{fi)C^{X)V^p,^p (54) 
S{X)Cpif,) = fifi - X)CMSiX) + 9if^ - X)C(,{X)S{fi) (55) 

taf3{X)Cy{fl) = /(A - n)Cp{fl)taaiX)V^p^f3y + g{X - fl)Cf3{X)ta^{fl) (56) 

/(A) = ^, 9W = -^,, (57) 

where V is the permutation operator for C"^ C"^. It is important to notice that the 
Latin and greek indices do not mix in these relations. 

One then applies the transfer matrix on the state |Ai, Xp) and with the help of the 
above relations commutes it through the C^s. The contributions from S and are 



treated exactly as in [|13|. The contributions from Taa (a 7^ 1) vanish for p < M while 
|Ai,...,Ap) is an eigenvector when p = M, without any constraint. To see this we need 
the commutation relations between the Cartan subalgebra generators and the creation 
operators C. Using relations (pi|), ( PB| ) and (^Tj), one easily derives 

[Hl\C,iX)] = -C,{X) (58) 
[H{\Cp{X)] = 0, \/aeA-{l} (59) 

Relations (^) imply that the eigenvector Ansatz has no a-states in it (a 7^ 1), while 



relation ( ^8|) implies that Cai{Xi)...Ca^{Xp) ||1) has only i?-states in it when p = M. This 
means 

C«,(Ai)...C,^(Ap) 111) = for p>M (60) 

Taa(A)C,,(Ai)...Ca^(Ap) 111) = 0, 0<p<M-l and a^l (61) 

Ta,(A)C,,(Ai)...C,,,(AM) 111) = (xsin(A))^^C,,(Ai)...C,,,(AM)||l), a ^ 1 (62) 

One then finds the corresponding eigenvalues of r(A) 

j^iin,,n,),M)^X) = (cos(A))^' f[ f (Xj - A) + (x"^ sin(A))^' ( n /(A - Aj) I A^"''''^ 

+ (rii-l)(xsin(A))npAf (63) 

Here A*^"^.?) jg eigenvalue of the unit-shift operator r^^^'^\ for a chain of p sites and 
^2 possible states at each site; it is constructed out of the permutation operator V on 
C"^ 0C'^^. The coefficients F°'^''"'°'^ are such F is an eigenvector of r'^"^'^) for the eigenvalue 
A("'2,p)_ Finally the Bethe Ansatz equations are just 

(-ir^f^^^^)'' = A^"-) ,J = 1,...,P (64) 
V sm(Aj) ; 
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The operator 

r("2-p) can be written as a product of disjoint permutation cycles. One also 
has l^r^"'^.?)^ = I _ Tj^g eigenvalues A*^"^'^) are then, at most, p^^ roots of unity and are 
highly degenerate. The dimensions of the cycles and their multiplicities will depend on 
both p and n2- 

One can perform the above diagonalization procedure over the pseudo- vacuum \\a) 
{a 7^ 1). The set of eigenvalues is exactly the same as the one found above and the 
eigenvectors have the same structure but form a completely distinct set, at least for 
< p < M. This is easily inferred from the action of the Cartan generators. One can 
also start with any of the vectors P G B and obtain yet other sets of eigenvectors. 
The superscript n2 is replaced by ni in (^), A^'^^'^^ and t^"-^'P\ These features reflect a 
large degeneracy of the spectrum. 

Finally, for n > 2, the eigenvectors |Ai,...,Ap) are generically not eigenvectors of all 
the magnetic field operators H^* |13[. Because the spectrum is degenerate this is not 
in contradiction with the fact that the magnetic operators commute with the conserved 
quantities. 

We have thus diagonalized the transfer matrix. The nested Bethe Ansatz truncates 
and the spectrum is trivial in the above sense. 

8 Conclusion 

A new algebra defining bosonic integrable 'free-fermions' representations has been derived. 
These models were shown to have a highly degenerate and 'simple' spectrum and to possess 
local symmetries. Another distinguishing feature of these models is the possibility to 
couple any two of them in an integrable way. The algebraic structure at the root of this 
'fusion' has been put in a simple algebraic setting, thereby unifying and simplifying all 
the derivations. New representations were found, generalizing the known XX models [|l^ 



and the Hubbard-like models ill 



Finding a multiparametric deformation of su{n) interpolating between the XXZ mod- 
els and the models considered here would be interesting |T3[. This would provide a 
quantum group structure which would shed a new light on the models at hand. 

Integrable bosonic Bariev chains and their multicopy generalizations have been found 



| 29| , 0. These chains, in their fermionic formulation describe correlated hoppings of 
electrons on a chain. The bosonic chains are obtained by coupling, in a yet different way, 
two or more su{2) XX models. It is therefore likely that one can use the free-fermions 
models to build new, more general Bariev chains. Martins reported some progress in this 
direction E^. The algebra found here is bound to describe the algebraic structure at 



work for the Bariev chains. 
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